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ABSTRACT 

In order to study the distribution of gas and ionizing radiation around rj Car and their impli- 
cations for its likely companion star, we have examined high-excitation emission lines of [Ne III], 
[Fe III], etc., in spectra obtained with the HST/STIS instrument during 1998-2004. Our principal 
results, some of them unexpected, are: (1) The high-excitation fluxes varied systematically and 
non-trivially throughout r\ Car's 5.5-year spectroscopic cycle. Instead of rising to a plateau after 
the 1998 "event," they changed continuously with a maximum in mid-cycle. (2) At one significant 
location a brief, strong secondary maximum occurred just before the 2003.5 spectroscopic event. 
(3) These emission lines are strongly concentrated at the "Weigelt knots" several hundred AU 
northwest of the star. With less certainty, [Ne III] appears to be somewhat more concentrated 
than [Fe III]. (4) A faster, blue-shifted component of each feature appears concentrated near the 
star and elongated perpendicular to the system's bipolar axis. This structure may be related to 
the equatorial outflow and/or to dense material known to exist along our line of sight to the star. 
(5) Using the photoionization program Cloudy, we estimated the range of parameters for the hot 
secondary star that would give satisfactory high-excitation line ratios in the ejecta. T c g rj 39000 
K and L ~ 4 x 10 5 L Q , for example, would be satisfactory. The allowed region in parameter 
space is wider (and mostly less luminous) than some previous authors suggested. 
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( Gaviolalll953l lAller fc Dunham|[l96l iThackera 



Eta Carinae is by far the most observable very 
massive star (M > 100 Mq), the only supernova 
impostor that can be observed in great detail, 
and it has other extraordinary attributes (se e 



1989; 



Hamann et al.l 199 



w 
4: 



19671: IViotti et al 

IZethsonl ~ 2001 ). Satisfactory analyses of these 



various papers m iHumphrevs fc Stanekl l2005h . 



The star and its ejecta produce more than 
2000 identified emission lines, mostly of com- 
plex singly ionized species such as Fe + and Ni + 
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low-ionization features have not been feasible, 
because they depend on radiative excitation pro- 
cesses which are extremely intricate. A few higher- 
ionization spectral lines, however, offer valuable 
clues because they are less complicated: [Ne III], 
[Ar III], [Fe III], etc. They originate in slow- 
moving gas that was ejected from 77 Car 100 to 
200 years ago. For reasons explained below, they 
help to constrain the nature of a hot companion 
star which is thought to exist, but which has not 
been observed directly. They also give clues to the 
nature and distribution of ejecta within 3000 AU 
of the central starQ 



1 Throughout this paper, if we omit individual references 
for either a minor detail or a well-known generality about 
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In this paper we use Hubble Space Telescope 
(HST) data to explore the variability and spa- 
tial distribution of the high-ionization lines. Sev- 
eral notable results emerge, including maps of the 
emission, continuous evolution of the intensities 
through 77 Car's 5.54-year spectroscopic cycle, and 
a suggestive blue-shifted component of each line. 
Some of these facts contradict models that pre- 
viously seemed credible. We also estimate limits 
to the secondary star's temperature and luminos- 
ity, based on the high-ionization features via pho- 
toionization calculations. We find a considerably 
broader range of allowed va lues than those pro- 
posed bv lVerner et al. ( 20051 ). 



the Weigelt knots, a set of thre e condensations 



1.1. 



Background information and astro- 
physical motivations 



Several distinct types of spectra occur within 
1" of 77 Car. (1) The stellar wind produces broad 
emission and absorption lines, several hundred km 
s _1 wide. (2) Slow- moving ejecta at r > 0.15" 
(i.e., more than 300 AU from the star) emit fea- 
tures narrower than 50 km sT 1 . (3) The same 
ejecta contain dust which reflects the stellar-wind 
spectrum. (4) Faster ejecta also show both in- 
trinsic and reflected spectra, though these are rel- 
atively faint in the central region. The high- 
excitation features belong mainly to category 2, 
but their excitation mechanisms are quite differ- 
ent from the narrow low-ionization lines. The lat- 
ter are strongly influenced by radiative processes 
such as resonance absorption of a UV photon lead- 
ing to the emission of two or more lower-energy 
photons. The high-ionization lines, however, de- 
pend on the same processes that dominate in or- 
dinary nebulae - especially photoionization by a 
hot star. This star is most likely the companion 
object mentioned in §1.21 below. 

The high-ionization lines have not been spa- 
tially mapped prior to this paper. Most of them 
are narrow and thus do not represent the 300- 
1000 km s _1 stellar wind. As a comparison, 
let us recall what is known about the spatial 
distribution of the low- excitation features. The 
first HST spectra of 77 Car showed that the 
narrow [Fe II] lines originate in the vicinity of 



n Car, see [Humphreys fc Staneld j2005f):lGull et~all j200lfl 



iMorse et al.l dl99Sh ; iDavidson fe Humphreys! [119971 ), "and 
numerous refs. cited therein. 



0.1" -0.3" northwest of the star (Davidson et al. 
19951 Il997t IWeigelt fc Ebersbergerlll986l) . These 



knots move outward rather slowly (v ~ 40 km 
s _1 ), pr obably near the equatoria l plane of the 
system (IHofmann fc Weigelt! 1 1988t I Weigelt et al. 
19951 iDavidson et al I Il997t IZethson et aIT Tl999: 



Dorland et all 120041 : ISmith et al.1 l2004h ). Later 



Space Telescope Imaging Spectrograph (HST/STIS) 
data confirmed that low-excitation lines col- 
lectively account for much o f the integrated 



brightness in that region (e.g. Gull et al. 1999; 



IJohansson et alJl2000tTZethson et alJll999l) . STIS 



could not resolve an individual Weigelt knot, but 
in cases where its slit crossed one of them, the 
spatial distribution of [Fe II] brightness peaked at 
the knot's position. Meanwhile many HST im- 
ages showed the individual knots at various wave- 
lengths, albeit not well resolved. Combining all 
these facts, most authors agree that the narrow 
low-excitation emission lines are largely concen- 
trated in the Weigelt knots. The high-ionization 
features, on the other hand, had little influence on 
the HST images and on Weigelt 's original speckle 
interferometryQ Given their very different ioniza- 
tion states and excitation mechanisms, it is not 
safe to assume that they coincide with the knots. 
They might represent diffuse lower-density gas 
surrounding the knots. This uncertainty inspired 
the main work reported below, §31 

1.2. Connections with 77 Car's 5.54-year 
cycle, and the secondary star 

Observers noticed long ago that these high- 
excitation lines occasionally weaken or even 
disappear and then gradually return to their 
norm a l strength after a few months ( Gaviol 



198 



1953t iThackeravl Il967t IWhitelock et all 
Zanella et al.l 1984 ). Zanella et al. interpreted 



2 ISmith et all (2000) asserted that HST/WFPC2 images 
with niters F631N and F658N showed [S III] and [N II] 
emission, respectively, in the Weigelt knots. However, (1) 
the F658N images were overexposed in the region of in- 
terest here, and also included an indeterminate amount of 
red-shifted Ha emission; and (2) there was no reason to as- 
sume that [S III] dominated in the F631N filter, which had 
been chosen to sample the continuum instead. STIS data 
show that in fact [S III] A6314 contributes less than 20% of 
the signal on the Weigelt knots. The WFPC2 F631N im- 
ages show primarily the fluctuating continuum rather than 
[S III] emission. 
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such "spectroscopic events" as shell ejections 
that temporarily quash the ionizing UV flux. 
Eventually enough events had been observed 
to show that the y recu r with a period of 5.54 



1994; 



years dDaminelil 19961 : Whitelock et al 
Damineli et all l2008bl ). The 2- 10 keV X-rays 



also d isa ppear on each o ccasion (jlshibashi et al 



Il999al |bl: ICorcoranl 120051 ). Many authors have 



discussed these phenomena, generally concluding 
that (1) there is a companion star in a highly ec- 
centric 5.54-year orbit; (2) the X-rays are formed 
by the two colliding stellar winds; and (3) a spec- 
troscopic event occurs near periastron due to a 
mass ejection and/or a disturbance in the wind 
and/or some sort of eclipse by the primary wind. 
The projected separation between stars is less than 
30 AU ~ 13 mas, unresolvable by HST, and the 
secondary star has not been directly detected (at 
least not in any u nambiguous way). For details 



SCO 



many r efs. i n Humphreys &: Stanekl ( 20051 ) . 



Gull et all (|200lh . and iMorse et al 



the phy sical considerations listed by 
(|2006bh . 



(1999), and 



Martin et al 



The simplest conventional theories do not pre- 
dict much variation of the high-excitation lines be- 
tween spectroscopic events. In those models the 
two stars are usually more than 10 AU apart and 
the densest regions of the primary wind are usually 
on the far side of the secondary, from our point of 
view. Therefore r\ Car's spectrum at average times 
has received comparatively little attention. In fJH 
though, we find interesting behavior between the 
1998.0 and 2003.5 events. 

The secondary star's quantitative parameters 
remain open to question. The opaque wind of the 
primary star probably h as a characteristic tem- 
perature below 20000 K (jrlillier et al.ll200ll ). too 
cool to account for the high-ionization lines. Its 
hypothetical companion, on the other hand, pre- 
sumably must be a massive star in order to cause 
an event, but not as massive as the primary and 
therefore less evolved; so we expect it to be hot. 
The X-ray spectrum impli es a 3000 km s _1 wind 
( Pittard fc Corcoranl 2002). and such a fast wind 
almost always implies T c s > 37000 K. 

The usefulness of a hot companion for exciting 
rj Car's high-ionization li nes has been recognize d 
for more than a decade (|Davidsonl H9971 fl999). 
This is an appealing idea because maybe we can 
deduce information about the secondary star from 



the high-excitation spectral features, via photoion- 
ization modeling. Since that object has not been 
observed directly, the 3000 km s^ 1 wind speed 
is almost ou r only other clue to its parameters. 
IVerner et al.l ( 20051 ) suggested spectral type and 
luminosity class 07.5 I, but it is not clear that 
such an object would fulfill the requirements. We 
explore this question in $5] below. 

1.3. Contents of this paper 

We restrict our analysis to a few well-understood 
emission lines that require ionizing photons with 
hv > 16 eV. Using spectra obtained with the 
HST/STIS thoughout an entire 5.54-year cycle, 
we examine the spatial distributions and temporal 
behavior of [Ne III] A3870, [Fe III] AA4659,4703, 
and the narrow component of He I A6680H Then 
we calculate photoionization models to match the 
relative intensities of [Ne III], [Ar III], He I, and 
other lines that require helium-ionizing photons 
(hv > 25 eV). 

In g2] we describe the HST/STIS data. The 
variability and distributions of the emission fea- 
tures are described in while 21 concerns an 
interesting and significant blue-shifted component 
seen in all of these lines. In fj5] we describe our 
photoionization modeling with implications for the 
secondary star, and |J6] concludes with a summary 
and brief discussion of our results. 

2. Observations, spectral analysis, and 
mapping technique 

Eta Car was observed with HST/STIS from 
1998.0 to 2004.2, covering an entire 5.54-year cy- 
cle. The first of these observations occurred on 
1998 January 1, a few weeks after the most crit- 
ical time in the 1998.0 spectroscopic event. The 
last of them were obtained 2004 March 6, about 8 
months after the 2003.5 event. Most of the obser- 
vations were concentrated in mid-2003. 

This instrument provided the highest spa- 
tial resolution ever attained in spectroscopy of 
rj Car. The STIS/CCD had three large-format 
(1024x1024 pixel) detectors with 0.05" square 
pixels, covering a nominal 52"x52" square field 
of view and operating from below 2000 A to 



3 Throughout this paper we quote vacuum wavelengths and 
heliocentric Doppler velocities. 
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about 10000 A. We used spectra observed with the 
52"x0.1" slit in combination with the G230MB, 
G430M, and G750M gratings which provided res- 
olutions R ~ 5000 to 10000. The observations 
include a variety of slit positions and orientations, 
with a concentration at position angle 332° where 
the star and Weigelt knots B and D all fell within 
the slit@ Altogether, a map of the slit locations 
resembles a quasi-random ensemble of intersecting 
lines (see the figures in ^3] below). 

HST/STIS/CCD data from the HST Treasury 
Project public archive^ include several improve- 
ments over the normal STScI pipeline and stan- 
dard CALSTIS reductions. Modeling and inter- 
polating the pixels in a special way improved the 
effective resolution along the cross-dispersion axis, 
comp ared to most o ther published STIS/CCD 
data (|Davidsonl l200(ih . In this process the data 
were rebinned so that one pixel in the reduced 
data corresponds to 0.5 original CCD pixel, about 
0.025". Initial bad/hot pixel removal, wavelength 
calibration and flux calibration matched or ex- 
ceeded the STScI pipeline and CALSTIS. We ex- 
tracted many local spectra centered at uniformly- 
spaced locations along the STIS slit; the interval 
between successive locations was 0.025" and each 
extraction was 0.1" wide along the slit. (Charac- 
teristics of the original data p recluded narrower 
extractions, see lDavidsonll2006l). At 77 Car's dis- 
tance D w 2300 pc dDavidson fe Humphreys! 19971 : 



tance L) sa pc ( Davidson tz riumpn 
Davidson etaH boOll ISmithl 120061 ). 



.1 corre- 



sponds to a projected size of about 230 AU. 

We examined the spatial distribution and the 
temporal behavior of [Ne III] A3870, [Fe III] A4659 
and A4703, and the narrow component of He I 
A6680. These were selected for minimal contam- 
ination by adjacent emission or absorption fea- 
tures. ([Ne III] A3870, for instance, is weaker than 
[Ne III] A3969 but the latter nearly coincides with 
bright hydrogen and helium features.) For the spa- 
tial maps we used data extending from November 
1998, 11 months after the 1998.0 event, to the 
beginning of June 2003, a few weeks before the 
2003.5 event (Table [2 we shall commment on the 



4 The star and the three main Weigelt knots are customarily 
labeled A, B, C, and D, with r e lative positions descr i bed b y 



Weieelt fc E bcrsbcrgcr ( 1986); Hofm ann fc Weieelj §988); 



Smith et aO l|2004t ). and marked in our Fig. [2] Their sepa- 



rations expand by about 1% per year. 
5 http:/ /etacar. umn.edu/ 



last dates below). We measured fluxes by inte- 
grating each line above the continuum, or in some 
cases above an underlying broad emission profile. 
Each measurement represented a spatial interval 
of about 0.1" along the slit, the extraction width 
mentioned above. 

The assembly of a spatial map for each emis- 
sion feature was a non-trivial exercise. Of course 
each STIS spectrogram sampled only the region 
covered by the slit, effectively a strip about 0.13" 
wide if we allow for HST's spatial resolution. On 
each observation date only one or two slit posi- 
tions had been used; but these varied in orienta- 
tion and some of them were offset from the star, 
so the whole ensemble sampled a useful fraction 
of the area within 0.5" of the star. Now consider 
any one of the measured emission lines. Subject 
to later verification, suppose that the feature's rel- 
ative spatial distribution did not vary much - i.e., 
we tentatively assume that maps at any two dif- 
ferent times closely resemble each other except for 
a normalization factor. (This statement does not 
apply to observations during an event, which we 
did not attempt to use.) Since the line's bright- 
ness varied through the spectroscopic cycle, we 
renormalized the values measured in each spec- 
trogram so that the feature in question always 
had an adjusted value of unity at the only place 
that was observed on every occasion: the cen- 
tral star. Fortunately some narrow-line gas ex- 
ists along our line of sight to the star, producing 
[Ne III], [Fe III], and He I lines superimposed on 
its spectrum. These reference measurements were 
less precise than those elsewhere, because the su- 
perimposed narrow emission features were weak 
compared to the underlying stellar spectrum. For 
each offset slit position, we adjusted the renor- 
malization factor to maximize the consistency at 
intersections with other slit positions. 

The above procedure obviously depends on 
some insecure assumptions, but we can apply con- 
sistency tests. Some of the slit positions passing 
through the star were close together or were used 
on more than one occasion, and each offset slit 
position intersected several of the others (see fig- 
ures in Ej3j> - Consequently a substantial number 
of points were independently re-observed at well- 
separated times. Note that each slit sample had 
only one adjustable parameter, its renormaliza- 
tion factor. Based on comparisons of the points 
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re-observed in separate spectrograms, the overall 
self-consistency of each map turned out to be quite 
good, and we found no local inconsistency worse 
than the expected measurement errors. Therefore 
our assumptions were sufficiently valid H Inciden- 
tally, possible time-delay effects cannot strongly 
influence these maps, since the light-travel-time 
for 0.5" is only about a week at the distance of 
r\ Car, and the cooling time in a He + region is 
less than a month for electron densities n e > 10 5 
cm~ 3 . 

In this paper, "phase" in the 5.54-year cycle 
is that used in the rj Car HST Treasury Pro- 
gram archive: P = 2023.0 days, * = MJD 
50814.0 = J1998.0 exactly. This zero point corre- 
sponds to phas e ~ 0.0 09 in a system proposed by 



Damineli et al 



(l2008bh . Zero points based on spe- 
cific phenomena in the spectroscopic events should 
be avoided for two reasons: The 1998.0-2003.5- 
2009.0 events differed in major respects, and, also, 
successive revisions of to chosen in this way have 
made comparisons among papers more difficult. 
All proposed orbit models are far too imprecise to 
be useful in this regard. The Treasury Program 
phase definition, on the other hand, has been ex- 
tant without alteration for a number of years and 
its period is consistent with observations. Con- 
sidering that the Treasury Program archive is the 
largest available set of data on n Car, it is the 
obvious reference standard. 

3. Behavior and spatial distribution of the 
emission 

The different photon energies required for 
[Ne III], [Fe III], and He I emission allow us 
to probe a range of physical parameters. He I 
A6680 is a recombination line formed in a He + re- 
gion, while [Ne III] and [Fe III] arise in Ne ++ and 
Fe ++ zones. Ionization potentials of H°, Fe + , He , 
Ar+, and Ne+ are 13.6, 16.2, 24.6, 27.6, and 41.0 
eV respectively. Thus, for reasons explained by 



6 These statements are true even for the data obtained a few 
weeks before the 2003.5 event. If more observations had 
been made in the mid-cycle period 1999—2002, we would 
have cautiously excluded the pre-event 2003 data; but for- 
tunately the latter turn out to be consistent. On 2003 May 
17 and June 2 the offset STIS slit intersected several earlier 
slit locations, see Table[T]and figures in jj3] After the renor- 
malization factor had been applied, line strengths matched 
reasonably well at those intersections. 



Osterbrock fe Ferland ( 2006f ). Fe ++ occurs mainly 
in regions of H + and He while Ne ++ and Ar ++ 
usually coexist with He + . 

3.1. [Ne III] emission line at 3870 A 

The top panel of Fig. [T] shows the time- 
dependent flux of the narrow [Ne III] A3870 emis- 
sion line that was seen superimposed on STIS 
spectra of the star (at resolution ~ 0.1") through 
the spectroscopic cycle 1998-2003. After a min- 
imum lasting several months, the line intensity 
increased slowly until mid-cycle and then grad- 
ually decreased again. The intensity curve did 
not form a plateau. Then, a few months be- 
fore the 2003.5 event, the line intensity rapidly 
increased again before falling off toward its min- 
imum. This late-cycle rise resembled the behav- 
ior of the He II emission and the observed X-ray 
light c urve 



(|lshibashi et all Il999bt iMartin et al 



2006 of ): perhaps the high-excitation emission at 
that time was somehow related to the X-rays 
or to the colliding-wind s hocked zones (see §6 



Humphreys et al. (|2008h for reasons why this 



might occur). 
Some, but not 



all, of th i s beha vior pattern 



was noted by iDamineli et alj ([2008a) for [Ar III] 
A7138. They measured equivalent widths (not 
fluxes) in ground-based spectra including the star 
plus ejecta out to r ~ 1". Their Fig. 1 differs from 
ours in two interesting respects: (1) The declines 
after mid-cycle do not match. At phase 0.8, for in- 
stance, the high-excitation lines had already fallen 
by 50% in the STIS results but less than 20% in 
the ground-based data. (2) Damineli et al. did 
not observe a brief, dramatic flux increase around 
phase 0.9. Most likely these discrepancies resulted 
from the very different spatial resolutions. The 
ground-based spectra include an amorphous un- 
resolved mixture of emission regions, whereas our 
Fig. prefers to a well-defined 0.1" locale which is 
particularly significant for a reason noted below. 

A simple model would have predicted an ex- 
tended plateau in Fig. [T] The [Ne III] flux shown 
there is superimposed on the star and has a neg- 
ative Doppler velocity of about —40 km s _1 ; so it 
must represent slow-moving gas located approxi- 
mately between us and the star. If it was ejected 
50-200 years ago, this gas has now moved several 
hundred AU from the star. If the orbit orienta- 
tion is within the range favored by most authors 
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(jlshibashil I2OO1I lOkazaki et all l2008f) . then dur- 
ing most of the 5.5-year cycle the initial part of 
a path from the hot secondary star to the Ne ++ 
in question (i.e., toward us) should pass through 
the low-density secondary wi nd facing away from 
the pr imary star; see Fig. 2 in lPittard fc Corcoran 
(j2002h . We omit quantitative details here, but pre- 
vious data would have allowed a scenario in which 
the ionizing photons are not seriously depleted 
along such a path until they reach the observed 
Ne ++ region. In that case the [Ne III] brightness 
superimposed on the star would change little dur- 
ing most of the cycle. Fig.Q]contradicts this simple 
model. 

One can easily imagine qualitative explana- 
tions for the observed behavior after seeing it, but 
choosing the right one is harder. Most likely the 
ionizing photon path mentioned above does inter- 
sect a substan tial and varying part of the dense 
prima ry wind. iGull et al. (2009) and lParkin et al 
(120091) . :'or instance, have discussed possible "con- 



ventional" models for the flow of dense gas out- 
ward from the binary system, allowing a substan- 
tial column density in the relevant sense. Soker 
(priv. comm.), on the other hand, remarks that 
Fig.[T]is qualitatively consistent with a very differ- 
ent picture, in which the secondary star is usually 
on the far rather than the nea r side of the primary 
(|Kashi fc Sokeri I2008L l2009bl ld\ A discussion of 
these issues is beyond the scope of this paper, but 
we note two points: (1) Every proposed model de- 
pends on a number of assumptions which are not 
easy to verify, and (2) the information in Fig. [1] 
is essential for any realistic view of the problem. 
Instead of an ill-defined average of ejecta around 
the star, this figure represents emission between 
us and the star at high spatial resolution. 

Fig. [2] shows a spatial map of the narrow 
[Ne III] A3870 feature in slit positions observed 
with STIS. In order to compensate for time varia- 
tions, we renormalized the measurements so that 
the value is always unity at the position of the star 
(see ^2]) . Most of the flux originates in the region of 
the Weigelt knots, at distances 0.1"-0.3" from the 
star. The positions of the rather amorphous knots 
BCD are marked in the figure. STIS could only 
marginally resolve them individually, but the peak 
obviously corresponds to Weigelt knot C. In our 
map we can locate this peak at r w 0.21", position 
angle w 300°. This corresponds well with mea- 



surements of knot C in HS T images wh i ch ga ve 
r w 0.22" in 1999-2003 (e.g. ISmith etHl (|2004) V 

Higher flux in [Ne III] A3870 was also observed 
at the locations of Weigelt knots B and D at po- 
sition angle ~ 335°. Their peak intensities are 
roughly half that of knot C. However, it is not pos- 
sible to distinguish between B and D based only 
on Fig. H 

Given the spatial-resolution limits of all HST 
data, the identification of [Ne III] with the Weigelt 
knots in Fig. [2] is practically as good as the evi- 
dence for their association with [Fe II] reviewed 
in i jl.ll above. Evidently the narrow-line [Ne III] 
emission does not chiefly represent a diffuse halo 
enveloping all the Weigelt knots, as seemed possi- 
ble before. 

3.2. [Fe III] emission lines at 4659 A and 
4703 A 

We examined both [Fe III] A4659 and A4703 
in order to compare them for mutual consis- 
tency. Their theoretica l intensity rat io is 1.85 
(|Nahar fc Pradhanl 1 19961: iQuinetl l'l996h . Figs. [3] 
and S] show maps of the narrow [Fe III] A4659 and 
A4703 fluxes respectively. In the same manner 
as for [Ne III], we rescaled the measured [Fe III] 
fluxes so the net value at the location of the star 
was always unity. Before rescaling, our measured 
flux of [Fe III] A4659 at that position was about 
1.9 times as high as for [Fe III] A4703 in each ob- 
servation, close to the theoretical value. The two 
[Fe III] maps are in excellent agreement with each 
other. [Fe III] emission originates in a slightly 
larger region than [Ne III]. This is not very sur- 
prising, given the difference in ionization poten- 
tials; [Ne III] occurs where helium is singly ionized 
but [Fe III] originates in zones where hydrogen is 
ionized but helium is not. 

Incidentally, just as this paper was completed 
we noticed a fainter additional [Fe III] condensa- 
tion about 0.48" NNE of the star, in data obtained 
in June 2009 as part of a STIS re-commissioning 
exerciseQ We mention this knot because it has 
a positive Doppler velocity, about +27 km s" 1 
rather than the values around —40 km s _1 seen 
in BCD. Its location was not sampled by the slit 
positions used in 1998-2004. 



7 HST program 11506: K. S. Noll, B. E. Woodgate, C. R. 
Proffitt, k, T. R. Gull. 
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In addition to [Ne III], Fig.[TJshows the tempo- 
ral variation of [Fe III] AA4659,4703 on the posi- 
tion of the star. They display the same qualitative 
behavior as [Ne III] A3870. 

This paper is not directly concerned with low- 
excitation features, but we must acknowledge that 
some of the narrow [Fe II] lines also showed tempo- 
rary minima at phases around 0.8 in Fig. [TJ qual- 
itatively like [Ne III] and [Fe III]. They all showed 
a peak at phase 0.9 and a brief minimum during 
the 2003 event. The [Fe II] lines varied differently 
from each other, but none of them showed a strong 
growth around phase 0.3 as [Ne III] did. 

3.3. He I emission line at 6680 A 

He I A6680 is a recombination line formed in 
He+ (|Osterbrock fe Ferlandll2006h . In a spectrum 
of rj Car, this feature usually consists of a broad 
component formed in the stellar wind, plus the 
narrow component from slow-moving ejecta. With 
the spectral resolution of STIS (~ 40 km s _1 ), the 
underlying broad emission makes the narrow He I 
line far more difficult to measure than [Ne III] and 
[Fe III]. Consequently it is not a clear enough in- 
dicator of the time dependence to include in Fig. 
[TJ Our He I A6680 map (Fig. [5j) is thus of lower 
quality than the others. There is no evident dis- 
agreement compared to [Ne III] . 

4. Blueshifted components of the high- 
excitation lines 

In addition to the narrow lines discussed above, 
[Fe III] and [Ne III] also show separate, broader 
components extending between Doppler velocities 
—250 km s^ 1 and —400 km s _1 . These also ex- 
ist for other high-excitation lines such as [Ar III], 
[S III], [Si III], etc. In each case the blue-shifted 
component is much wider than the narrow lines 
but much narrower than the normal stellar wind 
emission features, (rj Car's brightest Balmer lines 
usu ally have FWHM sa 4 00 to 500 km s _1 ; see, 
e.g., Davidson et al.ll2005h 

iThackeravi <|l953L Il967l ) recognized [Fe III] and 
[S II] components as emission blue-shif t ed b 
about -300 km s _1 . lAller fe Dunham! (|l96 
listed these and similar blue-shifted Fe II and 



ticularly including [Ne III]. They concluded that 
the blue-shifted emission arose at larger distances 
from the star; but the Weigelt knots and po- 
lar/equatorial morphology were not known at that 
time. Comparing the 19 92 spectroscopic event 
to the 1995 midcycle state, Damineli et al. ( 1998f ) 
confirm ed the Zanel l a et al . statements about vari- 
ability. iGull et all (|2009h recently showed STIS 
spectrograms (not tracings) of some of these fea- 
tures, but did not produce spatial maps like those 
presented below. 

Quantitative measurement of the blue-shifted 
[Ne III] A3870 is difficult due to nearby Si II A3864 
and Cr II A3867 (|Zethsonll200ll ). We therefore con- 
centrated on the blue-shifted component of [Fe III] 
A4659, centered at about 4653 A. To verify that 
it is [Fe III] emission, we compared the profiles of 
[Fe III] A4659 and A4703. Fig. [6] shows that the 
blue-shifted feature is undoubtedly a component 
of each [Fe III] line. This figure also illustrates 
the complexity of the line profiles, tracing various 
ejection features. There is probably some material 
with intermediate Doppler shifts. 

The temporal evolution of the blue-shifted com- 
ponent follows the spectroscopic cycle in roughly 
the same way as the narrow lines. The normal 
and the shifted components appear and disappear 
together. Fig. [7] compares the [Fe III] A4659 emis- 
sion line profile in spectra of the star at differ- 
ent phases in the cycle. The solid curve refers to 
a phase of 0.59, close to the time when the nar- 
row high-excitation emission was strongest. The 
dashed curve is the profile shortly after the 2003.5 
event, at phase « 1.0. The narrow component 
at 4659 A disappeared completely, and the broad 
blue component was diminished. 

Similar to fj3j in Figs. [8] and [9] we map the flux 
of the blue-shifted component of [Fe III] A4659. 
Most of it originates in the inner 0.1" region, less 
than 250 AU from the star. The distribution is 
not as sharp as a stellar point-source would be. 
The emission region is detectably elongated to the 
northeast and southwest, perpendicular to the axis 
of the bipolar Homunculus eject a- nebula H 



[Fe II] lines. IZanella et alJ |l984j) noted that some 
of them disappeared along with the ordinary nar- 
row high-excitation lines during an event, par- 



3 lGull et al] 1120091 ) mentioned "diffuse arcs" WSW and ENE 
of the star, but they meant velocity vs. position corre- 
lations in spectrograms, not spatial arcs. "WSW" and 
"ENE" referred to the position angle of the STIS slit when 
those obervations were made. Although they used direc- 
tion names such as NE and ENE, those authors did not 
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Considering the proximity of the blue compo- 
nent to the position of the star, its size, shape, 
and orientation, we suspect that Fig. [9] traces the 
near side of a latitude-dependent structure in the 
outer stellar wind. Its observed blueshift of about 
—300 km s _1 is smaller than the wind speed seen 
at th e bottoms of the P Cyg absorption fea tures 
(e.g JSmith et al.ll2003l : [Davidson et al.ll2005h . but 
these two quantities sample the gas in very differ- 
ent ways. Electron densities cannot greatly exceed 
10 7 cm -3 in the relevant gas, because higher val- 
ues would collisionally de-excite [Fe III], [Ne III], 
and [Ar III] too strongly (§5.3 below). Assuming 
that i] Car's mass-loss rate has the conventional 
value of ~ 10~ 3 M© yr" 1 and that its equatorial 
wind is l ess dense than its p olar wind by a factor of 

may fall below 10 7 
0.05". At 



200 



n., 



order 3 (jSmith et al 

~ 3 at equatorial radii r > 100 AU 



cm 



higher latitudes the corresponding radius would 
be larger. Since these size scales are much larger 
than the 5.54-year orbit, very likely the dense, bi- 
conical inner polar wind (r < 100 AU) shields 
the outer, lower-density high-latitude zones from 
the secondary star's ionizing radiation. In other 
words, at high latitudes the ionizing radiation is 
probably absorbed by gas that is too dense to pro- 
duce strong forbidden lines. Thus one expects the 
high-excitation forbidden lines to be concentrated 
toward latitudes less than, say, 50° or so - result- 
ing in an oblate emission region. If this picture 
is valid, evidently red-shifted emission from the 
far side cannot penetrate through the configura- 
tion. Local inhomogeneiti es or "clumping" d o not 
alter the basic reasoning. Gull et al. ( 2009h have 
recently discussed relevant time-dependent struc- 
tures at smaller size scales. Their models depend 
on a number of assumptions and do not make use 
of a polar wind; but in general they have little ef- 
fect on our comments above. Arguably the most 
significant well-defined observational clue is the 
elongation shown in Fig. [9]- a new result. 

On the other hand, we can imagine two other 
locations for emission farther from the star. First, 
the —300 km s _1 Doppler shift is well-matched 
to the region where our line of sight to the star 
intersects the inner parts of the bipolar Homuncu- 
lus ejecta-nebula and/or the "Little Homuncu- 



lus" (jlshibashi et alj 120031 ). The main objec- 
tion is that one might then expect the emission 
zone to have a projected size of 0.5" or larger. 
The other possibility involves a remarkable ob- 
served fact: our line of sight to the star has 
about two magnitudes mor e extinction than our 
view of the Weigelt knots ([D avidso n et al. 1995; 
Hamann et"aD Il999t Iffillier et al.1 l200lh . Dust 
grains cannot exist much closer to 77 Car than 
r ~ 150 AU (iDavidson fe Humphrevs1ll997[) . but 
apparently some very localized dusty material lies 
between us and the star. Can its far side or inner 
side be the region where the blue-shifted emission 
lines originate? Both the size scale in Fig. [TJ] and 
the Doppler velocity seem reasonable. If this is 
the explanation, then the blue-shifted line com- 
ponents must be intrinsically much brighter than 
the observable narrow lines but are obscured by 
the dust. (The same statement also applies to 
emission in the outer wind. This model and the 
outer-wind idea mentioned earlier are not neces- 
sarily distinct from each other.) 

5. Photoionization modeling and the Na- 
ture of the Secondary Star 

As mentioned in §1, the emission lines discussed 
here probably indicate photoionization by a hot 
star. Consider, for instance, [Ne III] A3870. This 
is a familiar feature in moderately high-excitation 
photoionized nebulae, where Ne ++ typically ex- 
ists in a He + zone maintain ed by ionizing photons 
with 25 eV < hv < 54 eV (|Osterbrock fe Ferland 
20061) . Within that zone, [Ne III] emission is col- 
lisionally excited by thermal electrons. Suppose 
we attempt to construct a non-photoionized model 
for the ejecta of 77 Car. Since Ne + has an ion- 
ization potential of 41 eV, the creation of suffi- 
cient Ne ++ by thermal collisional ionization would 
require some mechanism - e.g., low-speed shock 
waves - that can heat a substantial amount of 
material to T ~ 20000-30000 K. This is difficult 
to achieve because a single low-speed shock car- 
ries insufficient energy; and no such mechanism 
has been suggested by other data[f] Moreover, as 



describe the quantitative elongation shown in our Fig. [9] 
Their discussion focussed instead on velocity-position-time 
effects that one can see directly in the spectrograms. 



9 Shock fronts with v ~ 30 km s — 1 can produce the required 
temperature but the gas then cools rapidly. With densities 
n e ~ 10 5 to 10 8 cm -3 relevant in this paper, insufficient 
material would be hot unless it is shocked repeatedly, with 
a time interval less than a few months if n e ~ 10 5 cm -3 or 
a few hours if n e ~ 10 s cm -3 . A similar argument applies 
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Zanella et al.l (1984) emphasized, [Ne III] behaves 
very differently from the low-excitation lines dur- 
ing each spectroscopic event. Thus, as a very prob- 
able working hypothesis, we assume that [Ne III] 
X3870 signals ordinary, relatively straightforward 
quasi-nebular photoionization. So far as we know, 
no author has argued against this assumption in 
the case of r\ Car. If correct then it is most likely 
true for [Ar III] and [Fe III] as well. The low- 
excitation features, by contrast, are strongly influ- 
enced by local absorption in various UV permitted 
lines - "radiative pumping" - which depends crit- 
ically (and differently for each transition) on lo- 
cal velocity dispersions and density gradients. In 
other words the high-excitation lines are theoreti- 
cally more tractable than Fe II, [Fe II], and other 
low-excitation features. 

But where do the Ne + -ionizing photons above 
41 eV come from? Conventional models for the 



prim ary stellar wind are too cool (jHillier et al 



20011 ). The secondary star, however, is expected to 
have a high effective temperature. Its likely abil- 
ity to create He + and Ne ++ in the ejecta was rec- 
ognized as soon as its existence became probable 



(|Davidsonlll997l [T999). Since Fe++, He+, Ar++, 



and Ne ++ depend on stellar photon fluxes just 
above 16, 25, 27, and 41 eV, the relative strengths 
of high-excitation lines may indicate the energy 
distribution of the star's ionizing UV photons, and 
thereby its temperature. Since the glare of the pri- 
mary star makes this object impossible to study 
directly, practically the only other quantitative in- 
formation on it comes from the X-ray spectrum 
and from evolutionary constraints. Therefore, 
here we employ quasi-nebular photoionization cal- 
culations to model the relative high-excitation line 
strengths. Differences from earl ier calculations de- 
scribed by IVerner et al. (|2005l) will be noted in 



5.1. Procedure and assumptions 

We used versi on 08.00 of the pho toionization 
program Cloudy (jFerland et al. 1998) to simulate 
conditions observed at Weigelt knot C. Since the 
primary star can safely be neglected in the high- 



to heating by MHD waves. Strong shocks faster than 300 
km s — 1 are known to exist near r\ Car but they ionize neon 
far beyond Ne++. Other considerations and quantitative 
details are too lengthy to discuss here. 



ionization zones, we assumed a simple configura- 
tion with one stellar source and a uniform-density 
cloud at distance r « 10 16 cm « 700 AU which is 
reasonable for knot C. We chose a covering fac- 
tor of 0.05 and a filling factor of unity for the 
gas. We neglected possible UV extinction by dust 
within the He + region for three reasons: (1) local 
reddening of the [Fe II] emi ssion is known to be 
small in the Weigelt knots ( Hamann et al.l fl999t 



Davidson et al. 1995) while the high-ionization 
zones should have smaller column densities; (2) 
dust grains assoc iated with r\ Car tend to be 



abnormally large ( Davidson fe Humphreys! 1997h 
and therefore relatively ineffective at hu > 20 eV; 
and (3) not enough information is available to real- 
istically include the effects of local dust on the ion- 
izing radiation. Altogether our assumptions repre- 
sent an idealized view, but fortunately the line ra- 
tios used here depend only weakly on density and 
on geometrical details unless n e exceeds 10 7 cm -3 . 
A higher-density regime appears unlikely because 
collisional de-excitation would then suppress the 
forbidden lines too much. If n e < 10 7 cm -3 , the 
results depend mainly on the star's effective tem- 
perature T e ff and the ionization parameter E/h, see 
g3] below. 

We assumed that the He + zone is ionization- 
limited, i.e., that it absorbs nearly all incident 
photons above 25 eV. This assumption is based on 
three observed clues: (1) lower-ionization regions 
clearly exist in and around the Weigelt knots, (2) 
almost no He I, [Ne III], and [Fe III] emission is 
seen at radii outside the locations of the knots, 
and (3) ratios of the He I, [Ne III], and [Fe III] 
lines do not vary much. 

Studies of i] Car's outer ejecta show a pecu- 
liar chemical composition due to the CNO cycle. 
The gas is helium- and ni trogen-rich while car- 
bon and oxygen ar e scarce ( Davidson et al.lfl986t 
Dufour et al.l Il999h F°l We considered two sets of 



abundances, "compositions A and B" which dif- 
fer by a factor of two for elements heavier than 
helium (Table [2]). The mass fraction of C+N+O 
(mostly N) is about 0.63% in composition A and 
1.3% in B. The former value resembles some crude 



observational estimates for r\ Car (jDavidson et al 



10 The stellar wind spectrum appears consistent with the 
outer ejecta but does not indicate abundances as well 
jHillier et al.ll200lfl . 
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19861) . while the latter is somewhat larger than So- 
lar or Galactic material. The distinction between 
A and B turns out to affect the results only weakly, 
as explained in ^5.31 Note that carbon and oxygen 
are too scarce to play an appreciable role. 

We chose Weigelt knot C for this analysis be- 
cause the high-excitation emission is strongest 
and best-defined there, see Fig. [2 Knot C was 
observed well with HST/STIS on 2003 Febru- 
ary 13 (phase = 0.92 in the 5.54-year cycle) and 
2003 May 17 (phase = 0.97). We used data from 
2003 February 13, obtained with STIS gratings 
G230MB, G430M and G750M and with the slit 
passing through the star and knot C along posi- 
tion angle 303°. On that date the high-excitation 
lines were still strong, whereas by May 17 they had 
declined seriously as the 2003.5 event approached. 
We measured well-isolated narrow high-excitation 
lines as target values for the photoionization mod- 
els. The resulting equivalent widths and apparent 
intensities are listed in columns 4 and 5 of Table[3] 
An error of ±30% in an individual line strength, 
or perhaps even worse, would have little effect on 
our conclusions - partly because the photoioniza- 
tion models are very sensitive to the star's T e g, 
and partly because the He/Ne/Si/Ar abundance 
ratios are more uncertain than the observed line 
ratios. 

We estimated the reddening from the [Ar III] 
A3110/A5193 and [Ar III] A7138/A7753 flux ra- 
tios. The intrinsic value of the former ratio, for 
instance, is reliably known because the A3110 and 
A5193 lines share the same upper level; and the 
case of A7138/A7753 is similar. Other reddening 
indicators are less trustworthy because they use 
more than one ion species or differing upper levels 
or other model-dependent factors. Since no physi- 
cal model is available for r\ Car's anomalously large 
grains or for instrumental effects (see below), we 
adopted a standard approximation for wavelength- 
dependent extinction: A\ w a + b/X. For small 
amounts of reddening this choice of mathematical 
form is not critical, provided that coefficient b is 
adjusted to give the right average slope from violet 
to red wavelengths. Fitting it to the [Ar III] mea- 
surements, we found surprisingly little reddening, 
Eb-v = Ab — Ay « 0.2 magnitude. 

This is only half the a mount that the [ Fe II] 
lines seem to ind icate ( Davidson et al.l Il995l : 
Hamann et al.l[l~999h . but instrumental effects re- 



lated to high spatial resolution probably account 
for most of the discrepancy. The normal wave- 
length dependence of HST's spatial resolution 
caused the slit throughput to decrease toward 
longer wavelengths for a localized so urce, an d 



other effects also occurred in STIS data (Davidson 



2006). Since these complications had smooth 



wavelength dependences, they are implicitly in- 
cluded in the effective or apparent reddening de- 
duced from the [Ar III] line ratios. In other words, 
the Eb-v value mentioned above was really the 
true interstellar and circumstellar value minus a 
correction for instrumental effects. Meanwhile 
it is also possible that the [Fe II] method gives 
an overestimate of £b-v- No matter which effect 
dominates, the corrected relative line strengths are 
automatically valid to sufficient accuracy because 
the [Ar III] comparison method is based on known 
intrinsic ratios. (As mentioned earlier, "sufficient 
accuracy" in this context would be ±30% or al- 
lowably even worse. Most of the values in Table [3] 
are expected to be better than this.) 

Thus, based on the [Ar III] lines, we adopted 
the following correction for net reddening: 



-^corr = 0.105 exp 



/ 3740 A 



*obs ■ 



Here the constant factor, which has no effect on 
the line ratios used in our photoionization anal- 
ysis, is adjusted to give [Ne III] A3870 a cor- 
rected value of 1.00. As a result the corrected line 
strengths I corr are intrinsic values relative to this 
line. They are listed in the last column of Table 
® 

5.2. Stellar model atmospheres 

We explored a multidimensional grid of pho- 
toionization models, varying the effective stel- 
lar temperature, ionization parameter, and gas 
density. Our goal was to constrain these prop- 
erties by comparing the calculated intensity ra- 
tios of high-excitation emission lines to the ob- 
served ones. We tried four theoretical stel- 
lar atmosphere grids available in Cloudy. The 
Atlas models are LTE, plane-parallel, hydro- 
static atmospheres wi th turbulent veloc i ty dis - 
tribution 2 km s" 1 (|Castelli fe Kuruczl |2004 . 
The CoStar O-type models are non-LTE, line- 
blanketed model atmospheres, i nclud ing stel- 
lar winds ( Schaerer fe de Koter 1997f ). The 
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Tlusty models are non-LTE, line- blanketed, plane- 
paral lel, hydrostatic O-star SEDs (|Lanz fc Hubenv 
2003). The WM-basic O-star grids represent non- 
LTE, line-blanketed, wind-blanketed hot stars 
(|Pauldrach et al.ll200ll ). 

Unfortunately these four types of theoretical 
models disagree with each other in their UV spec- 
tral energy distributions. Fig. ITOl shows the con- 
tinuum of an O-type main sequence star with 
T cff = 40000 K and L = 10 5 L Q according to each 
model. Their differences are worst at high photon 
energies, particularly above 40 eV. Table[4]lists the 
total luminosities and photon rates that can ion- 
ize H°, He , and Ne + in a comparable set of the 
four model types. The CoStar models do not in- 
clude line opacity and therefore overestimate the 
far-UV flux because photon line-absorption and 
subsequent re-emission at longer wavelengths is 
not taken into account. The lower flux in the He 
continuum (hv > 25 eV) of the WM-basic model 
compared to the Atlas models is probably due to 
non-LTE effects producing deeper line cores in the 
blocking lines. The WM-basic code uses a con- 
sistent treatment o f line blocking and blanketing 
( Smith et al. 2002). After experimenting with all 
four types, we adopted this type of model. 

5.3. Resulting constraints on the sec- 
ondary star 

The Weig elt knots encompass a wide range 
of densities ( Hamann et alj Il999h . We there- 



fore varied the hydrogen density from 10 5 to 
10 7 cm~ 3 , and for each density we varied the 
ionization parameter £/h between 10 -2 and 10 2 . 
In Cloudy, Uh is defined as the dimensionless 
quantity QH/47rr 2 riHC, where Qh is the rate of 
hydrogen-ionizing photons (hv > 13.6 eV) emit- 
ted by the source star, r is the distance from that 
source to the illuminated face of the cloud, nn is 
the hydrogen density, and c is the speed of light 
For our purpose, of course, we care most about the 
photon supply above 25 eV rather than 13.6 eV. 
Cloudy produces ionization-limited plane-parallel 
models, obviously a crude approximation to the 
true geometry. 



11 The ionization parameter determines the sharpness of 
each ionization front, the coexistence of differing ioniza- 
tion stages, the ionized col umn density, and other ph ysical 



Based on these calculations, Tableland Fig.fTTI 
show the effective temperatures of WM-basic stel- 
lar atmospheres that give the observed strengths 
of Si III] A1892, [Ar III] A7138, and He I A6680 
relative to [Ne III] A3870. ([Fe III] and [S III] are 
satisfactory in the models favored below.) For ion- 
ization parameters logt/n > +1 the required stel- 
lar temperature is roughly constant, while smaller 
values of logC/ H require progressively higher T e g. 
Since these results involve a subtle blend of un- 
certainties, the shaded regions in Fig. [IT] repre- 
sent factor-of-4 ranges for each line ratio. The 
best match to the observed data set occurs with 
log/7 H < -1 and T cS ~ 40000 K. Here are some 
relevant considerations: 



1. 



It is not surprising that chemical composi- 
tions A and B ( <j5. 1[) give similar results. 
Roughly speaking, the hydrogen and he- 
lium recombination lines indicate numbers 
of ionizing photons that have been ab- 
sorbed, while the heavier-element emission 
lines account for much of the cooling and 
thus represent th e total energy in the ab- 
sorbed photons (iDayidson fc Netzer 1979t 
Osterbrock fc Ferlandl 120061) . Therefore, if 



attributes, as explained by iDavidson &: N ctzcr (197! 



we alter the overall abundance of heavy ele- 
ments relative to H+Hc, the equilibrium gas 
temperature automatically adjusts so that 
the ratio of heavy-element emission lines to 
hydrogen and helium lines does not change 
much. This ratio depends chiefly on the 
slope of the ionizing source's spectral energy 
distribution. 

2. For \ogUu > —1, He I A6680 indicates gen- 
erally lower values of T e g than the other 
lines do. Very likely this clue is evidence for 
logC/jr < —1; but we also note that models 
which are density-limited, or limited by in- 
ternal dust, or convex, can produce stronger 
helium lines relative to [Ne III] . 

3. Since relative abundances of individual 
heavy elements are quite uncertain in rj Car, 
and the reddening correction is uncertain for 
the UV Si III] line, the systematic difference 
between [Ar III] and Si III] in Fig. QT] is not 
very alarming. 

4. Although the high-ionization emission lines 
are easier to model than low-ionization fea- 
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tures, some special processes may occur. For 
example Si III] A1892 may be enhanced by 
reso nance-favored two-pho ton ionization of 
Si + ( Johansson et al. [20061 ). Such effects are 
unusual, however, and unlikely to alter our 
basic conclusions. (If the Si III] effect is 
strong, it produces a discrepancy between 
He I and Si III] in photoionization models.) 

5. Models with n e > 2 x 10 7 cm~ 3 produce 
insufficient [Ne III] and [Ar III] emission 
compared to He I and Si III], regardless 
of the stellar temperature. This is because 
the characteristic electron densities for col- 
lisional de-excitation of [Ne III] and [Ar III] 
are roughly 10 7 and 5 x 10 6 cm -3 respec- 
tively. 

6. Models with n e << 10 6 cm -3 are ruled out 
by geometrical considerations. For reasons 
mentioned in H5.1[ the He + region cannot 
be much smaller than the ionization-limited 
thickness. This characteristic linear size x 
is proportional to U-n/n e . In Tableland 
Fig. [11] we have indicated models that are 
unsuitable because x > 250 AU, i.e., an 
ionization-limited He + zone would be larger 
than the region of the Weigelt knots. This 
criterion excludes all our models with n e = 
10 5 cm -3 , and those which have n e = 10 6 
cm~ 3 and logf7 H > -0.65. 

7. A truly realistic model would include a range 
or distribution of gas densities, but we don't 
have enough observables to do this. At 
present one can say only that the "repre- 
sentative density" in a simplified model is 
very likely about half the maximum density 
which exists in the real, non-uniform He + 
gas, in an order-of-magnitude sense. (This 
statement can be wrong if the configuration 
is unexpectedly complex.) 

8. For a given gas density n e and stellar tem- 
perature T c ff, one can estimate the lumi- 
nosity L that would produce the assumed 
value of C/r in Weigelt knot C. As shown in 
the last column of Table [5] some choices of 
(n e , Uh) lead to absurdly small or large val- 
ues of L. Most important, the secondary star 
cannot exceed about 10 6 Lq since its pres- 
ence is not evident in the UV spectra ob- 



tained with HST/STIS. (The primary star 
has L w 4 x 10 6 i Q .) 

The H-R diagram in Fig. [12] summarizes our 
conclusions about the secondary star, assisted by 
a few more clues. Here two broad curves mark 
the results of our photoionization calculations with 
n e = 10 7 cm~ 3 and 10 6 cm -3 ; the upper end of 
the latter curve is limited by criterion 6 above. If 
we assume that the two stars have the same age, 
then evolutionary circumstances further constrain 
the parameters. The minimum age for the pri- 
mary star to have become helium- rich is 0.5 Myr 
(|lberJll999h . but it is most likely 2-3 Myr old based 
on its instabilities and probable associ ation with 
the cluster Trumpler 16 (|Walbor nlll995h . The life- 
time of any star above ~ 60 Mq is roughly 3 Myr. 
Therefore, in Fig. [12] we show isochrones for 0.5, 
2, and 3 Myr as well as some evolution ary tracks, 
all adapted from iMartins et al.l (J2005J) . We can 
base yet another constraint on the observed X- 
ray spectrum, which indicat es a secondary wind 
speed close to 3000 km s _1 (jPittard fe Corcoran 
2002). Practically all stars with winds that fast 
have T e g > 370 00 K, marked by a vert ical dotted 
line in Fig. [TJ (jKudritzki fc Pulsll2000l and refer- 
ences cited therein). 

Altogether, then, we expect the secondary star 
to lie somewhere within the shaded polygon in Fig. 
[711 with the following limits: (1) The right-hand 
boundary is the X-ray-implied minimum temper- 
ature. Even without this argument we would de- 
duce a similar limit from criterion 6 listed above 
for our photoionization models. (2) The upper 
limit is L f» 10 6 i©, criterion 8. (3) The upper- 
left boundary is based on our photoionization cri- 
terion 5, n e < 2 x 10 7 cm~ 3 . (4) The lower-left 
boundary is the 0.5 Myr isochrone, and perhaps 
this limit should be moved rightward in the di- 
agram to some age greater than 1 Myr. A star 
with M ZA ms ~ 40 to 60 Mq and T eS w 40000 K, 
for example, would satisfy all these requirements. 
The corresponding spect ral type would be 04 or 
05 ( Martins et al.l 120051 ). If, however, the sys- 
tem is more than 2 Myr old, then the secondary's 
zero-age mass was probably less than 50 Mq. (All 
statements related to isochrones obviously depend 
on the evolution models, though.) 

If the entire region around r\ Car were clearly 
visible, then the absolute brightnesses of the high- 
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excitation emission lines would indirectly indicate 
the secondary star's luminosity. In fact the local 
extinction is far too patchy, but we attempted a 
crude estimate based on the observable [Ne III] 
brightness. In our Cloudy calculations with rea- 
sonable input parameters, the ionizing star's lu- 
minosity was around 2300 times the luminosity of 
[Ne III] A3870. Measuring the total flux of this 
line in the Weigelt knots, and assuming reasonable 
factors for the extinction and the solid angle in- 
tercepted by the knots, we found logL/L© ~ 5.5. 
This estimate is very rough, but it is consistent 
with Fig. [12] and thereby suggests that the abso- 
lute fluxes are reasonable. 

IPittard fc Corcoran! (|2002f) assumed that the 
secondary star has a mass-loss rate close to 
10~ 5 Mq yr _1 , in order to obtain sufficient X-ray 
luminosity in their colliding-wind models. Such a 
high rate would be very unusual, and not entirely 
consistent with our discussion above. Possibly this 
is an argument in favor of a more luminous sec- 
ondary star; but on the other hand the mass-loss 
estimate is not very robust. Only a small fraction 
of the wind's kinetic energy is converted to ob- 
servable X-rays, via some complicated efficiency 
factors. The secondary wind speed, however, al- 
most directly determines the average temperature 
seen in the 2-10 keV X-ray spectrum. This is why 
we consider the 3000 km s _1 speed estimate to 
be relatively more trustworthy than the mass-loss 
rate. 

One can imagine models that are seriously af- 
fected by changes in the secondary star. For in- 
stance, conceivably it was originally more massive 
but h as now become a Wolf-Rayet object. ISoker 
(l2007h has proposed another scenario, wherein the 
secondary star accreted a large amount of mass 
during the great eruption 160-170 years ago, and 
has not yet retur ned to its normal therm al equi- 
librium. (See also lKashi k. Sokerl (|2009al U Apart 
from obvious complications and a lack of substan- 
tive evidence for them, models of these types are 
clearly not the simplest possibilities. They have 
multiple adjustable parameters, and thus are be- 
yond the scope of this paper. (Note, incidentally, 
that our photoionization calculations do not favor 
an extremely hot star with with T e g > 45000 K.) 



5.4. Comparison with previous calcula- 
tions 



Verner et alj (|2005l ) reported an earlier set of 



photoionization calculations very much like ours, 
with the same goal; but they deduced apprecia- 
bly different parameters for the secondary star. 
Like us, they used Cloudy with WM-basic atmo- 
spheres, uniform density in the ionized material, 
etc. They used earlier STIS data, but the high- 
excitation line ratios were similar to ours within 
the uncertainties. Verner et al. appear to have 
supposed that the high-ionization lines originate 
much closer to the star than the locations of the 
Weigelt knots, but this probably had little effect 
on their results. They described photoionization 
models with specific sets of parameters, but did 
not show systematic maps of parameter space as 
in our Tableland Figs. E] and HU 

There are three notable differences between 
their results and ours: The allowed region in Fig. 
1121 is much larger than their discussion seems to 
imply, their suggested upper limit for T e g is barely 
above our lower limit, and they assigned a max- 
imal luminosity to the secondary star. Indeed 
Verner et al. proposed that it has T e g w 37200 
K and logL/L© ~ 5.97, marked with a small "V" 
in our Fig. [T2 (They clas s ified it "07.5 I", but 
according to lMartins et alJ (|2005l ) an 05.5 super- 
giant would have that temperature.) Such an ob- 
ject is not excluded by the calculations, but in 
our view it has three disadvantages: It has the 
minimum temperature required for a 3000 km s _1 
wind, it requires the age of the system to be less 
than 2 Myr, and it has practically the largest al- 
lowable luminosity. So far as we know, there is 
no strong argument against lower luminosities. In 
summary, given the limited information currently 
available ( §5.31 above), the parameters suggested 
by Verner et al. are not the most suitable choices 
for assessing the nature of r\ Car's companion at 
this time. 

6. Summary 

Here we have reported the behavior and the 
spatial origins of rj Car's quasi-nebular high- 
excitation emission between spectroscopic events. 
Neither of these questions was previously explored 
with adequate spatial and temporal sampling. In- 
deed the 1998-2004 HST/STIS spectra appear 
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to constitute the only suitable existing data set. 
Thus our main results were previously vague or 
unknown, and some of them contradict models 
that would otherwise have seemed credible. 

For instance, after the 1998 spectroscopic event 
the high-excitation features did not simply recover 
to constant "normal" intensities, but instead var- 
ied systematically throughout the 5.5-year cycle 
(§523 GT5J). [Ne III] and [Fe III] evolved through 
four distinct stages (Fig. [TJ: 

1. As expected based on previous events, the 
weak or near-zero state persisted for several 
months. This was long enough for a compan- 
ion star to sweep through more than 180° of 
longitude around periastron, if its orbit ec- 
centricity is at least 0.8 as seems likely. 

2. Subsequent growth to near-maximum in- 
tensity was interestingly slow, delaying the 
maximum until more than two years after 
the event - unlike what a simple model 
would predict. 

3. Then, instead of leveling off, the intensity 
soon began to decline, following a sort of 
parabolic trajectory in time. 

4. A few months before the 2003.5 event, how- 
ever, it briskly rose to a second maximum 
just before declining again to near-zero. 

We directly measured the double-peaked cycle in 
gas along the line of sight to the star, but at least 
stages 1-3 seem broadly valid in our data on the 
Weigelt knots BCD a s well. As noted in §3.1| 
Damineli et a l. ( 2008a) recognized those stages for 
one [Ar III] line at much lower spatial resolution, 
but there is a quantitative disagreement in stage 
3 and they did not see stage 4. 

If this time-pattern was caused by some com- 
plex variation of the circumstellar extinction, then 
in order to include BCD it must extend out to 
r > 800 AU, which seems unlikely in terms of grain 
physics. Apparently, then, during most of the 
1998-2003 cycle the amounts of Ne+ 2 , Fe+ 2 , etc. 
in BCD were significantly less than their maxima. 
Despite our comments at the beginning of fj5j the 
secondary maximum might conceivably be related 
to the colliding-wind shock structure which was 
very strong in the months before a n event. The 
softest X-rays may play a role then (jMartin et al 



2006bl ). A detailed speculative/theoretical in 
terpretation of Fig. [T] would requir e calculations 
beyond the s cope of this paper. Martin et al 



(<2006aL l2009h lave described other phenomena 
observed midway between rj Car's spectroscopic 
events, emphasizing that its mid-cycle behavior 
deserves more attention. 

In §3 we also showed spatial maps of the [Ne III] 
and [Fe III] emission in the sub-arcsecond vicinity 
of i] Car. Informal examinations of [Ar III], [S III] 
and [Si III] gave similar impressions. Evidently 
high- excitation emission is strongly concentrated 
in the Weigelt knots BCD located at a distances 
0.15"-0.35" northwest of the star, rather than in a 
surrounding halo or in the region between the star 
and the knots as seemed possible before. Prior to 
this work most researchers would very likely have 
assumed that each knot produces high-excitation 
emission, but here we have given the first clear 
proof. Weigelt blob C has the highest flux. Maps 
of the various ion species all show basically the 
same picture, except that those with higher ioniza- 
tion potentials seem a little more compact. Unfor- 
tunately it will be difficult to improve these maps 
in the future, because the rising brightness of the 
central s tar is rapidly makin g the knots harder to 
observe (jMartin et al.ll2006d) . 

The high-excitation lines all have blue-shifted 
components in the velocity range —250 to —400 
km s . This fact was known before, but in §4 
we have shown that ( 1 ) it originates mostly within 
0.15" of the central star, and (2) this region is elon- 
gated along the equatorial (not polar) axis of the 
rj Car system. In §4 we proposed several plausi- 
ble interpretations involving the three-dimensional 
distribution of stellar ejecta. These alternatives 
are not mutually exclusive. 

Assuming that a hot companion star is respon- 
sible for the narrow high excitation emission lines, 
in Sj5]we found constraints on its parameters based 
on photoionization calculations. Our results are 
consistent with evolutionary considerations and, 
independently, with the colliding-wind X-ray tem- 
perature. Fig. [12] shows that the allowed region 
in parameter space is larger than some previous 
authors suggested. For example, an 04-06 gi- 
ant with L ~ 4 x 10 5 L Q , T off » 39000 K, and 
-Mzams ~ 40-50 Mq would fall near the center of 
the allowed range. Our calculations did not ex- 
clude moderately higher masses and luminosities, 
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but neither did they favor them. 
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Fig. 11.— Contour plots of line intensity ratios [Ar III] A7138/[Ne III] A3870 (top), He I A6680/[Ne III] A3870 
(middle), and Si III] A1892/[Ne III] A3870 (bottom) in our photoionization models with nn= 10 6 cm~ 3 . 
Shaded areas correspond to the range of line intensity ratio between 0.5 and 2 times the observed value. Left 
and right columns refer to chemical compositions A and B respectively. Models above the dashed curve are 
spatially too extended, violating criterion 6 in the text. 
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Fig. 1.— Strengths of the narrow [Ne III] A3870 
and [Fe III] AA4659,4703 emission lines measured 
in HST/STIS spectra of r) Car through its spectro- 
scopic cycle. These represent gas along the line of 
sight to the star, seen with 0.1" spatial resolution 
and not including the Weigelt knots. Features in 
the knots vary in a similar way. 
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Fig. 2.— Spatial map of the narrow [Ne III] A3870 
line. Position A marks the central star while B, 
C, and D are the Weigelt knots. Flux values are 
normalized so that this emission line's net flux is 
unity at the star's position. 
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Fig. 3. — Spatial map of the narrow [Fe III] A4659 
line, with the same format as Fig. [2] 
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Fig. 4.— Similar to Fig. HH1 but showing [Fe III 
A4703. This map is based on the same spectro- 
grams as Fig. [3] but the measurements were inde- 
pendent. 
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Fig. 6. — Profiles of the two [Fe III] emission lines 
that we measured at the position of the star. Ar- 
rows mark the broad negative-velocity emission 
and vertical ticks show the heliocentric zero-shift 
wavelengths. 




0.15 -0.15 -0.3 -0.45 
x in arcsec 



Fig. 5. — Map of the the narrow He I A6680 emis- 
sion line, shown in the same way as Figs. [2HH 
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Fig. 7. — Comparison of the [Fe III] A4659 emis- 
sion components at two different phases in the 
5.54-year cycle. The solid line represents a phase 
of 0.59 near high-excitation maximum, while the 
dashed line shows phase 1.00 during the 2003.5 
event. 
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Table 1 

Times of STIS observations used in the emission maps 



Date (UT) 


MJD 


Slit PA a 


JV B m b 


1998 Nov 25 c 


51142 


227° 


1 


1999 Feb 22 


51231 


332° 


1 


2000 Mar 21 


51624 


332° 


1 


2001 Apr 18 


52017 


22° 


1 


2001 Oct 01 c 


52183 


165° 


1 


2002 Jan 20 


52294 


278° 


1 


2002 Jul 05 


52460 


69° 


2 


2003 Feb 13 


52683 


303° 


1 


2003 Mar 29 


52727 


332° 


1 


2003 May 05 


52764 


27° 


1 


2003 May 17 


52776 


38° 


2 


2003 Jun 02 


52792 


62° 


2 



a PA ± 180° gives same spatial coverage. 

The star was observed on every occasion, 
but sometimes there was also a second slit po- 
sition offset about 0.25 NW of the star. 

c There were no [Fc III] obs. on this date. 



Table 2 
Chemical compositions 



logio N(element)/N(H) 
Comp. A a Comp. B b 



He 


-0.70 


-0.70 


C 


-5.00 


-4.70 


N 


-3.10 


-2.80 


O 


-5.00 


-4.70 


No 


-4.00 


-3.70 


Si 


-4.46 


-4.16 


S 


-4.74 


-4.44 


Ar 


-5.60 


-5.30 


Fe 


-4.55 


-4.25 



a Default solar composition as used in 
Cloudy, except He, C, N, O. 

b CNO fraction somewhat larger than 
solar. 
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Table 3 

Selected high-excitation emission lines in Weigelt blob C 



Spectrum 


Avac 




I.P.'s a 


EW 


T b 

Jobs 


7 corl c 




(A) 




(eV) 


(A) 






[Ne III] 


3869.85 


40. 


.96-63.45 


10.24 


3.63 


1.00 


Si III] 


1892.03 


10 


34-33.49 


3.84 


2.78 


2.10 


[Ar HI] 


3110.08 


27 


63-40.74 


0.42 


0.17 


0.06 


[Ar III] 


5193.26 


27 


63-40.74 


0.39 


0.11 


0.024 


[Ar III] 


7137.76 


27 


63-40.74 


6.36 


2.24 


0.40 


[Ar III] 


7753.24 


27 


63-40.74 


1.93 


0.59 


0.10 


He I 


4027.33 


24 


59-54. 42 d 


2.42 


0.69 


0.18 


He I 


6680.00 


24 


59-54. 42 d 


5.44 


2.54 


0.47 


He I 


7067.20 


24 


59-54. 42 d 


16.19 


10.1 


1.80 


[S III] 


6313.81 


23 


33-34.83 


7.72 


2.39 


0.45 


[S III] 


9071.11 


23 


33-34.83 


4.06 


1.46 


0.23 


[S III] 


9533.23 


23 


33-34.83 


5.49 


2.49 


0.39 


[Fe III] 


4659.35 


16 


.18-30.65 


4.27 


1.97 


0.46 


[Fe III] 


4702.85 


10 


.18-30.65 


3.21 


0.88 


0.21 



References. — Atomic data from fhttp:/ /physics. nist. gov /j Phys- 
RcfData/ASD/. 

a Rclevant ionization potentials. For instance, Nc + and Ne"^~*" 
have ionization potentials 40.96 and 63.45 cV. 

b In units of 10 — 12 erg cm" 2 s _1 within the 0.1" sampled area. 

c Dc-rcddcncd intrinsic strength relative to [Nc III] A3870. 

d Rccombination spectrum created in He"^ zone. 



Table 4 

Comparison between different stellar atmosphere models' 1 



Model 


L 

(erg/s- 1 ) 


log L/Lq 


log Qh 
(photons s _1 ) 


log Qhc 
(photons s _1 ) 


log QnbIII 
(photons s~ 1 ) 


Atlas 


3.67c+38 


4.98 


48.72 


48.27 


46.09 


CoStar 


4.25e+38 


5.04 


48.90 


48.70 


48.32 


Tlusty 


3.27o+38 


4.93 


48.76 


48.40 


46.33 


WM-basic 


3.99e+38 


5.02 


48.82 


48.51 


47.04 



a T off = 40000 K, n H = 10 6 cm" 3 and log Un = -1 
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Table 5 

wm-basic model predictions of t eff for secondary star and composition a 



Density 
(cm ) 


log Uk 


log Qh 
(photons s~ 1 ) 


T eff (Ar 7138) a 
( K ) 


T eff (Ho 6680)" 
( K ) 


T aff (Si 1892)" 
( K ) 


T 
(K) 


AT 

( K ) 


log L/Lq 


10 5 


-2.14 


46.44 


48700 


40300 


37600 


42200° 


5800 


2.84 


10 5 


-2 


46.58 


48700 


38900 


37200 


41600 b '° 


6200 


2.99 


10 5 


-1 


47.58 


42300 


35900 


35100 


37800 b,c 


3900 


4.08 


10 5 





48.58 


39100 


34200 


33900 


35700 b 


2900 


5.16 


10 5 


1 


49.58 


37700 


33300 


33500 


34800 b 


2500 


6.22 


10 5 


2 


50.58 


36800 


33000 


33100 


34300 b ' c 


2200 


7.31 


10 6 


-2 


47.58 


43400 


37600 


41100 


40700° 


2900 


4.00 


10 6 


-1 


48.58 


38100 


33900 


36700 


36500 


2300 


5.13 


10 6 


-0.67 


48.91 


37200 


33600 


37200 


36000 


2100 


5.48 


10 6 





49.58 


36600 


32800 


36600 


35300 b 


2200 


6.18 


10 6 


1 


50.58 


36200 


32600 


36300 


35000 b ' c 


2100 


7.19 


10 6 


2 


51.58 


36600 


32800 


35900 


35100 b ' c 


2000 


8.19 


10 7 


-2 


48.58 


39400 


38000 


51900 


43100 


7700 


4.97 


10 7 


-1 


49.58 


36100 


34900 


43700 


38200 


4800 


6.06 


10 7 





50.58 


36000 


34500 


38700 


36400 c 


2100 


7.13 


10 7 


0.24 


50.82 


35700 


34500 


38400 


36200° 


2000 


7.38 


10 7 


1 


51.58 


35500 


34900 


38100 


36200 b ' c 


1700 


8.14 


10 7 * 


2 


52.58 


35600 


34900 


37800 


36100 b ' c 


1500 


9.14 



a Rclativc to [Ne III] A3870 A. 

b Modcls with He II region larger than 4 X 10 15 cm 
c Modcls with log L/Lq smaller then 5 or larger than 6. 
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Fig. 8. — Spatial map of the relative flux in the 
blue-shifted component of [Fe III] A4659. The flux 
is normalized to unity at the position of the star. 
When measured carefully the bright region is elon- 
gated along the NE-SW direction, see Fig. [9l 
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Fig. 9. — Expanded view of intensity contours 
of the blue-shifted [Fe III] A4659, using the same 
measurements as Fig. [5] Squares (■) and triangles 
(A) indicate values of 0.5 and 0.25 respectively, 
relative to the value at the star's location. 
The contours are significantly elongated NE-SW, 
i.e., toward the upper left and lower right. 
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Fig. 10. — Comparison of the stellar continuum in 
Atlas, CoStar, Tlusty and WM-basic atmosphere 
models with T cff = 40000 K. Tlusty and WM-basic 
models were smoothed. 
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Fig. 12. — Likely positions for rj Car's secondary 
star in the H-R diagram. Solid curves show 
isochrones for 0.5, 2, and 3 Myr, while dashed 
curves show evolution track s for initial masses 1 5 
to 85 Mq, all adapted from lMartins et al.l (|2005r ). 
The two broad lines show L vs. T e g correlations 



in our photoionization models with tih = 10 7 and 
10 6 cm" 3 . A small 'V indicates a model suggested 



bv lVerner et al. (2005). The shaded polygon is the 
region allowed by various considerations (see text) . 
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